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摘 要 : 准确 估计 输 水 条 件 下 河岸 地 下 水 埋 深 的 动态 变化 ,可 以 量化 生态 输 水 量 与 地 下 水 至 深 的 


上 的 合理 性 和 适用 性 。 然 后 根据 3 个 断 


具有 重要 的 科学 意义 。 结 合 塔 里 木 河 下 游 20 a4 
的 拟 二 维 地 下 水 模型 ,对 输 水 条 件 下 塔里木 河 下 游 上 、 中 、 下 段 3 个 断 间 
麻 ) 的 地 下 水 埋 深 变化 进行 20a 长 期 模拟 。 通 过 率 定期 和 后 11 a(2010 一 2020 年 ) 的 地 下 水 埋 深 模 
拟 结果 与 站 点 数据 比较 ,发 现 两 者 较 一 致 , 证 明 该 模型 在 塔里木 河 下 游 河岸 断 本 
20 a 的 模拟 结果 分 析 输 水 条 件 下 地 下 水 埋 深 和 土壤 水 的 
长 期 变化 及 其 对 生态 输 水 的 响应 。 结 果 表 明 : 经 过 20a 的 生态 输 水 , 英 苏 阿拉 干 和 依 干 不 及 麻 3 
个 断面 上 的 地 下 水 位 和 土壤 湿度 都 有 明显 的 上 升 ,地 下 水 位 至 深 从 输 水 前 的 8 四 左右 抬升 到 输 水 
后 的 4mm 左 右 ,土壤 温度 从 最 初 的 0.20 上 升 到 0.35 以 上 ,特别 是 自 2009 年 以 来 , 随 年 输 水 量 增 加 ， 
地 下 水 位 和 土壤 湿度 增加 幅度 明显 。 生 态 输 水 与 地 下 水 的 年 际 变 化 有 一 定 的 滞后 性 ,由 于 土壤 湿 
度 和 地 下 水 位 表现 为 正 相 关 关 系 , 这 使 得 土壤 湿度 对 输 水 量 也 有 灌 后 性 的 特点 。 相 比 于 河水 流 
量 ,地 下 水 水 平 传导 率 的 取 值 对 断面 地 下 水 埋 深 变 


响应 关系 ,并 由 此 估计 自然 河道 所 需 输 水 量 及 持续 时 间 , 这 对 于 干旱 区 水 资源 管理 的 可 持续 发 展 
E 态 输 水 监测 数据 ,用 发 展 的 包含 地 下 水 和 土壤 水 


( 英 苏 、 阿 拉 干 和 依 于 不 及 


地 下 水 长 期 模拟 


比 起 着 更 重要 的 作用 。 男 外 , 输 水 量 与 地 下 水 


的 年 际 变化 表明 塔里木 河 下 游 河 岸 要 想 获得 持续 的 生态 效益 ,需要 对 河道 提供 间歇 性 的 生态 输 水 。 
关 键 词 : 拟 二 维 地 下 水 模型 ; 生态 输 水 ; 地 下 水 长 期 模拟 ; 塔里木 河 下 游 
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塔里木 河 是 中 国 最 长 的 内 陆 河 ,流域 位 于 新 疆 
南部 。 塔 里 木 河 流域 在 历史 上 形成 了 天 然 绿洲 ,是 
阻挡 风沙 侵袭 .保护 人 类 生存 环境 的 重要 生态 屏 
障 中 。 然 而 , 自 1950 年 起 , 随 着 大 规模 农业 和 人 工 
绿洲 等 人 为 措施 的 实施 ,使 得 水 资源 被 过 度 开发 利 
用 ,导致 塔里木 河源 流 和 干流 水 资源 分 配 不 均 ,至 
使 塔里木 河 下 游 自 大 西海 子 水 库 至 台 特 玛 湖 321 km 
河道 断 流 , 河 两 岸 地 下 水 埋 深 下 降 , 植 被 大 片 死亡 ， 
土壤 沙化 更 加 严重 ,生态 环境 遭 到 严重 破坏 2” 。 为 
了 挽救 生态 环境 ,实现 可 持续 发 展 战略 , 自 2000 年 
起 ,水 利 部 同 新 疆 维 吾 尔 自治 区 人 民政 府 及 新 疆 生 
产 建设 兵团 等 部 门 对 塔里木 河 下 游 实施 了 间隔 性 
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生态 输 水 工程 ,特别 是 2010 年 以 来 ,区 域内 的 生态 
植被 得 到 了 有 效 的 恢复 “”。 了 解 和 掌握 塔里木 河 
沿岸 地 下 水 和 土壤 水 在 输 水 工程 实施 以 来 发 生 的 
变化 ,对 于 定量 评估 输 水 工程 的 生态 效应 ,以 及 干 
早 区 生态 水 资源 管理 提供 重要 的 科学 依据 "1。 
目前 ,针对 输 水 条 件 下 河岸 地 下 水 变化 研究 的 
主要 手段 包括 站 点 监测 和 模型 模拟 。 陈 亚 宁 等 
在 塔里木 河 下 游 河道 的 上 、 中 和 下 段 的 9 个 断面 上 ， 
根据 离 河岸 距离 的 远近 建立 了 多 个 地 下 水 埋 深 监 
测 井 , 用 来 分 析 地 下 水 对 河道 生态 输 水 的 响应 。 随 
后 ,许多 研究 "根据 这 9 个 断面 上 监测 井 的 水 位 
数据 ,分 析 了 不 同 生 态 输 水 时 期 下 ,断面 上 地 下 水 
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埋 深 的 变化 情况 。 由 于 监测 井 在 断面 上 是 离散 分 
布 的 ,因此 这 些 数据 只 能 反映 了 该 井 的 水 位 变化 ， 
对 于 未 建立 监测 井 地 方 的 水 位 无 法 估计 ,即使 结合 
站 点 数据 ,建立 简单 的 线性 回归 模型 进行 估计 Al 
也 可 能 存在 着 较 明 显 的 误差 。 

地 下 水 模型 能 结合 有 限 的 站 点 数据 ,模拟 到 除 
了 监测 井 外 断面 上 其 他 点 的 地 下 水 埋 深 ,是 对 站 点 
监测 的 一 种 有 益 补 充 。 杨 鹏 年 等 所 提出 了 用 一 维 
地 下 水 非 稳定 流 模型 ,解析 模拟 输 水 条 件 下 河道 断 
面 上 地 下 水 埋 深 的 变化 。 随 后 刘 迁 迁 等 中 和 古 力 
米 热 . 哈 那 提 等 " 虽 对 非 稳定 流 模型 进行 了 不 同 程 
度 的 改进 以 提高 模拟 效果 ,但 归根 结 底 还 是 通过 对 
模型 做 各 种 假设 把 非 线 性 模型 变 成 线性 模型 求 得 
解析 解 , 这 些 假设 带 来 的 误差 会 比 用 数值 解 的 误差 
BK, 与 之 不 同 的 是 , 狄 振 华 等 1 对 地 下 水 非 稳 定 
流 模型 进行 了 离散 的 数值 解 , 并 在 谢 正 辉 等 “的 土 
二 水 运 边界 理论 的 基础 上 ,发 展 了 土壤 水 和 地 下 水 
相互 作用 的 数值 模型 ,并 在 英 苏 断面 上 得 到 了 较 好 
的 地 下 水 模拟 。 然 而 ,该 研究 只 模拟 了 输 水 期 内 英 
苏 断 面 上 地 下 水 水 位 的 变化 ,并 没有 对 其 进行 长 时 
间 的 地 下 水 模拟 ,而 且 对 塔里木 河 下 游 其 它 河 段 上 
的 地 下 水 水 位 也 未 进行 模拟 比较 。 为 此 ,本 研究 结 
合 塔里木 河 下 游 20a 的 生态 输 水 监测 数据 ,对 塔 里 
木 河 下 游 的 上 .中 和 下 段 各 选取 一 个 断面 ,分别 进 
行 20a 的 地 下 水 埋 深 长 期 模拟 ,定量 评估 输 水 条 件 
下 河岸 地 下 水 和 土壤 水 的 变化 。 


1 研究 区 数据 与 方法 
1.1 研究 区 概况 


研究 区 位 于 塔里木 河 下 游 大 西海 子 水 库 至 台 
特 玛 湖 的 河道 (图 1) ,河道 全 长 321 km ,东部 为 库 鲁 


大 西海 子 水 库 
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注 :C3~C6 为 上 段 英 苏 断面 上 的 4 个 监测 井 ;G3~G6 为 中 段 阿 拉 干 断面 上 
的 4 个 监测 井 ;H2~H5 为 下 段 依 干 不 及 麻 断 面 上 的 4 个 监测 井 。 
图 1 塔里木 河 下 游 河道 及 9 个 断面 示意 图 


Fig.1 Schematic diagram of the lower reaches of 


Tarim River and its nine sections 


繁 改道 的 过 程 中 沉积 ,形成 了 约 350 m 的 第 四 系 松 
散 沉积 层 。 沉 积 物 以 黏 士 质 的 河 、 湖 相 沉 积 为 主 ; 
地 形 平坦 , 微 向 东 倾 斜 ,坡度 3%o 左 右 。 河 两 岸 受 风 
沙 的 影响 ,周边 的 土壤 多 以 粗 粉 沙 - 极 细 沙 - 细 沙 
为 主 3, 这 种 地 质 状况 下 地 下 水 与 土壤 水 的 作用 较 
为 强烈 。 由 于 受 河水 漫 溢 及 地 下 水 抬升 影响 ,植被 
主要 沿 河 呈 带 状 分 布 ,植被 类 型 自 上 而 下 总 体 呈 由 
高 覆盖 度 的 胡杨 林 、 树 柳 灌 从 . 盐 生 草 带 向 中 低 覆 
HE A aK PERRE AA Ty (ER at EY Pa REE IE A, 
Te EEE o 

自 1950 年 以 来 ,人 类 不 合理 利用 水 土 资源 使 得 
大 西海 子 水 库 以 下 河道 长 期 处 于 断 流 状 态 , 地 下 水 
位 下 降 到 8~12 m, 导 致 由 地 下 水 维系 的 天 然 植被 严 
重 退化 ,土壤 盐 渍 化 等 生态 环境 问题 日 益 突 出 。 为 
了 挽救 塔里木 河 下 游 的 生态 环境 ,水利 部 同 新 疆 维 


克 沙 漠 ,西部 为 塔克拉玛干 沙漠 ,海拔 高 程 在 
801.50~846.24 m 之 间 ,属于 冲积 平原 区 。 该 区 年 降水 
Ht 17.4~42.0 mm ,年 潜在 蒸发 量 高 达 2500~3000 mm, 
实际 年 蒸发 量 178.5 mm 左右 ,属于 典型 的 大 陆 性 干 
旱 气 候 , 生 态 环境 极为 脆弱 。 由 于 有 河流 和 地 下 水 
补给 , 沿 河 两 岸 生 长 了 大 片 的 荒漠 河岸 林 , 阻 止 了 
两 大 沙漠 的 合拢 ,因此 具有 重要 的 生态 意义 和 战略 

研究 区 地 处 罗布 泊 微 弱 抛 陷 区 ,构造 稳定 ,在 
漫长 的 地 质 时 期 内 ,河水 携带 的 大 量 泥 沙 在 河道 频 


在 尔 自治 区 人 民政 府 及 新 疆 生 产 建设 兵团 等 部 门 
从 2000 年 5 月 开始 在 大 西海 子 以 下 河 段 实施 生态 
输 水 工程 ,直到 2020 年 共计 21 次 生态 输 水 ,累积 输 
水 量 84.4x10 m。 为 了 评估 生态 输 水 条 件 下 地 下 水 
埋 深 和 土壤 水 含量 的 变化 ,本 研究 结合 20 a 生态 输 
水 监测 数据 ,选取 塔里木 河 下 游 上 段 英 苏 断面 .中 
段 阿拉 干 和 下 段 依 干 不 及 麻 3 个 河 段 断 面 , 分 别 进 
行 地 下 水 和 土壤 水 的 模拟 ,以 期 揭示 生态 输 水 对 河 
岸 地 下 水 和 土壤 水 变化 的 有 影响。 根据 断面 上 监测 
井 的 位 置 , 本 研究 选取 河道 单 侧 断 面 的 1000 m 范 围 
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区 域 进行 模拟 ,断面 和 主要 监测 井 的 示意 图 见 图 
1。 英 苏 断 面 的 4 个 监测 井 (C3、C4、C5、C6) 阿拉 干 
断面 的 4 个 监测 井 (G3、G4、G5 .G6) 和 依 干 不 及 麻 断 
面 的 4 个 监测 井 (H2、H3、H4、H5) 分 别 距离 河岸 
150 m,300 m,500 m,750 m> 
1.2 数据 来 源 及 预 处 理 

(1) APT Fi 20 a 流量 数据 来 自 于 新 疆 塔 
里 木 河 流域 管理 局 ;断面 河道 流量 及 地 下 水 埋 深 数 
据 来 源 于 中 国 科 学 院 新 疆 生 态 与 地 理 研究 所 ,我 们 
只 收集 了 前 5 次 输 水 (2000 一 2003 年 ) 条 件 下 的 部 分 
河道 流量 和 断面 监测 井 地 下 水 数据 ,其 中 地 下 水 数 
据 有 英 苏 汤面 C3、C4、C5、C6 监测 井 , 阿 拉 干 断面 
G3 G4., G5 监测 井 和 依 干 不 及 麻 断 面 H2、H3 监测 
井 。 男 外 ,2010 一 2020 年 的 监测 井 地 下 水 埋 深 数据 
也 来 自 中 国 科学 院 新 疆 生 态 与 地 理 研究 所 。 值 得 
注意 的 是 ,本 研究 所 采用 的 模型 是 个 依赖 初 边 值 的 
物理 模型 ,只 要 给 定 初始 地 下 水 埋 深 值 和 河道 的 每 
日 流量 作为 边界 值 ,就 能 开展 长 期 的 模拟 。 

(2) 地 面 高 程 数据 根据 断面 上 监测 并 的 地 面 高 
程 数 据 插值 得 到 。 它 被 用 来 计算 地 面 到 地 下 水 埋 
深 的 垂直 距离 ,用 于 垂直 土壤 水 运 边 方 程 的 计算 。 
土壤 参数 据 根据 研究 资料 和 查 表 得 到 。 马 晓 东 等 呈 
通过 地 质 调查 发 现 塔里木 河 下 游 的 土壤 质地 以 沙土 
和 沙壤土 为 主 ,它们 的 饱和 含水 量 在 0.45~0.50 mm”, 
我 们 根据 下 游 河 段 的 经 纬度 范围 ,查找 BATS 模式 
对 全 球 土壤 的 12 种 典型 分 类 表 ”2 ,最终 确定 该 研究 
区 为 第 6 类 沙 塘 土 : 饱 和 土壤 水 含量 (8) 为 0.48 mm”, 
饱和 基质 势 (y.) 为 -0.20 m, HE ET PS KR K) H 
0.54432 med ,土壤 质地 参数 (5) 为 6.00。 结 合 文献 资 
料 和 查 表 ,得 到 的 饱和 含水 量 基 本 一 致 ,可 以 认为 
所 确定 的 土壤 参数 是 合理 的 ,因此 ,它们 将 被 用 于 
研究 区 土壤 水 运动 方程 的 计算 。 
1.3 研究 方法 
13.1 拟 二 维 地 下 水 模型 ”本 人 研究 的 地 下 水 与 土壤 
水 相互 作用 模型 采用 狄 振 华 等 ”发 展 的 拟 二 维 地 
下 水 模型 。 该 模型 是 将 断面 上 完全 二 维 的 土壤 水 
与 地 下 水 相互 作用 问题 简化 为 以 水 平流 为 主 的 地 
下 水 运动 程 和 以 垂直 流 为 主 的 土壤 水 运动 方程 相 
耦合 的 拟 二 维 模型 。 具 体 包括 : 


oh _ ð oh 
n= A (Kad) OREN CA) ( 1 ) 
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式 中 : 球 是 饱和 给 水 度 mm5a; 玉 是 地 下 水 水 平 导 水 
率 (m*d);h(xt) 是 1: 时刻 离 河岸 x 米 处 的 地 下 水 海 拨 
高 度 (m); q- 0) 是 在 潜水 面 上 地 下 水 与 上 层 土 
壤 水 间 的 垂直 交换 通 量 (md');0 为 土壤 水 含量 也 称 
土壤 湿度 (mi m°); DOM KO HIEM ERK 
BCE (m? d) A EK KR (m+ d’), E MTE FB Clapp 和 
Hornberger?” 关系 式 与 饱和 土壤 水 含量 ( 6. ,mm )、 
垂直 饱和 导 水 率 (Km:d4) .饱和 基质 势 ( wy ,m) 和 
土壤 质地 参数 (4 ) 建 立 关 系 式 。 g(x,) 是 地 表 的 
入 渗 ( 蒸 发 ) 通 量 (md")。 本 研究 中 模型 的 地 表 通 量 
采用 多 年 平均 实际 蒸发 量 ,模型 采用 隐 式 差分 格式 
算法 进行 数值 计算 ,具体 计算 细节 见 狄 振 华 等 ”。 
1.3.2 SCE-UA 参数 率 定 方法 SCE-UA 方法 是 
Duan 等 天 提出 的 全 局 参数 自动 优化 方法 , 它 结合 了 
遗传 算法 和 单纯 性 法 的 优点 ,不 需要 对 模型 进行 直 
接 求 导 ,能 快速 找到 全 局 最 优 解 ,是 当今 水 文 模型 
参数 优化 的 主要 流行 算法 之 一 。 同 狄 振 华 等 2 的 
SEC-UA 用 法 一 样 , 本 研究 主要 优化 拟 二 维 模型 中 
的 地 下 水 水 平 导 水 率 (K.) 参 数 , 而 假定 非 饱 和 土壤 
中 垂直 饱和 导 水 率 (K.,) 不 变 , 其 与 BATS 全球 土 壤 
分 类 表 中 的 饱和 传导 率 值 一 致 。 

1.3.3 模拟 实验 设置 ”地 下 水 方程 (1) 的 左边 界 是 
河道 选取 变化 的 河水 水 位 作为 左边 界 值 ,右边 界 设 
为 等 值 透水 面 , 故 选 取 零 通 量 为 边界 。 土 壤 水 方程 
(2) 的 上 边界 是 地 表 , 选 取 多 年 平均 地 表 通 量 (降雨 
减 蒸发 减 径流 ) 为 边界 值 , 取 178.5 mm-a”, FUI 
选取 固定 饱和 土壤 水 含量 为 边界 值 , 取 0.48 mm >。 
地 下 水 埋 深 的 初始 值 为 模拟 第 1 天 的 观测 值 ,土壤 
水 的 初始 值 因 无 观测 , 故 设 潜水 面 上 方 1 m 均 为 
0.4,1 m 以 上 至 地 表 均 为 0.2。3 个 断面 的 模拟 期 都 
是 输 水 的 20 a(2001—2020 年 ), 但 具体 来 说 又 略 有 
不 同 。 英 苏 的 模拟 期 为 2000 年 11 月 16 日 一 2020 年 
12 月 31 日 ;阿拉 干 的 模拟 期 为 2002 年 8 月 15 日 一 
2020 年 12 月 31 日 ; 依 干 不 及 麻 的 模拟 期 为 2002 年 
9 月 23 日 一 2020 年 12 月 31 日 ,其 中 第 2 次 输 水 期 的 
2000 年 11 月 16 日 一 2001 年 2 月 4 日 的 81 d 为 英 苏 
断面 的 参数 率 定期 ,第 4 次 输 水 期 的 2002 年 8 月 15 
日 一 11 月 3 日 的 91 d 和 2002 年 9 月 23 日 一 11 月 20 
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日 的 59 d 分 别 为 阿拉 干 和 依 干 不 及 麻 的 参数 率 定 
期 ,其 他 年 份 为 无 观测 约束 的 模拟 期 。 率 定期 的 河 
道 日 流量 和 河水 水 位 由 观测 值 获得 ,其 它 模 拟 期 
(如 2003 一 2020 年 ) 的 河道 日 流量 数据 , 则 是 根据 率 
定期 的 总 输 水 量 与 河道 断面 累积 流量 的 比例 关系 ， 
结合 当年 的 生态 输 水 量 和 输 水 天 数 进行 插值 得 
到 。 其 他 模拟 期 的 河水 水 位 则 根据 率 定期 河流 量 
与 水 位 的 拟 合 关系 式 及 插值 的 日 流量 计算 得 到 。 
率 定期 第 一 天 的 地 下 水 埋 深 作 为 模拟 初 值 ,后 续 的 
20 a 地 下 水 埋 深 通过 模型 模拟 得 到 。 根 据 收集 的 3 
个 断面 监测 井 地 下 水 埋 深 时 长 ,选取 模拟 期 的 后 11 a 
(2010—2020 年 ) 为 模型 验证 期 。 根 据 监测 井 的 位 
置 ,本 研究 分 别 选取 河道 3 个 断面 单 侧 的 1000 m 范 
围 区 域 进行 模拟 ,模拟 的 空间 步 长 是 1 m, 时 间 步 长 
是 0.5 he 


2 结果 与 分 析 


2.1 拟 二 维 地 下 水 模型 参数 率 定 
结合 观测 的 地 下 水 埋 深 数据 ,用 SEC-UA 参数 


8340 (a1) 英 苏 -C3 井 
& 832 5 814 
A 5 = 
和 830 Š 812 


810 


8167 (b1) 阿拉 干 -G3 井 808 


优化 方法 ,分 别 对 各 自 断 面 上 的 拟 二 维 地 下 水 模型 
进行 地 下 水 水 平 导 水 率 ( 玉 ) 参 数 进行 率 定 ,最 终 英 
苏 、 阿 拉 干 和 依 干 不 及 麻 断 面 上 的 KK. 分别 被 率 定 为 
1.588 m.d .0.079 m-d?,0.171 md'。 为 了 检验 最 
优 率 定 参 数 的 合理 性 ,比较 分 析 了 断面 上 观测 和 包 
含 最 优 参数 值 模型 模拟 的 地 下 水 水 位 , 英 苏 断面 上 
4 个 监测 井 (C3、C4、C5、C6)、 阿 拉 干 断面 3 个 监测 井 
(G3 、G4、G5) 和 依 干 不 及 麻 断 面 2 个 监测 并 (H2、H3) 
在 率 定期 内 观测 和 模拟 地 下 水 水 位 变化 结果 如 图 2 
所 示 。 通 过 比较 发 现 , 在 经 过 参数 率 定 后 ,模拟 的 
地 下 水 位 高 程 与 其 观测 值 比较 接近 ,其 中 英 苏 4 个 监 
测 井 地 下 水 埋 深 模拟 的 平均 绝对 值 误 差 (MAE ) 为 
0.194 m, 阿 拉 王 3 个 监测 并 模拟 的 MAE 为 0.214 m 
和 依 干 不 及 麻 2 个 监测 井 模拟 的 MAE 为 0.189 mo 
同样 地 ,也 可 以 看 到 在 同一 断面 上 模拟 的 地 下 水 水 
位 离 河 越 远 而 变 得 越 低 。 这 些 结果 表明 , 拟 二 维 地 
下 水 模型 能 够 合理 模拟 塔里木 河 下 游 生态 输 水 条 
件 下 河岸 的 地 下 水 变化 ,从 而 证 明了 模型 在 该 区 域 
是 适用 的 ,为 下 步 的 长 期 输 水 条 件 河岸 地 下 水 模拟 
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图 2 3 个 断面 监测 井 在 率 定 


期 内 模拟 和 观测 的 地 下 水 水 位 变化 


Fig. 2 Variations of the simulated and observed groundwater table for the monitoring wells at the three sections 
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提供 合适 的 模型 平台 。 
2.2 断面 地 下 水 埋 深 的 长 期 模拟 

在 获得 3 个 断面 上 各 自 的 最 优 模型 参数 值 后 ， 
结合 大 西海 子 水 库 对 塔里木 河 下 游 的 20 a 生态 输 水 
量 数据 ,就 可 以 用 拟 二 维 地 下 水 模型 对 输 水 条 件 下 
河流 断面 上 的 地 下 水 埋 深 ,进行 20 a 长 期 模拟 。 根 
据 收集 到 的 2010 一 2020 年 断面 监测 井 地 下 水 位 ,我 
们 对 模拟 结果 和 观测 进行 了 比较 ,结果 如 图 3 所 
示 。 模 拟 值 是 断面 上 1000 m 内 的 地 下 水 埋 深 年 平 
均值 ,观测 值 是 断面 上 5 个 监测 井 的 年 均 地 下 水 埋 
深 值 ,其 中 每 个 断面 上 5 个 监测 井 选取 距离 河岸 
150 m,300 m,500 m,750 m, 1050 m 的 监测 井 。 最 
终 英 苏 断 面 上 的 地 下 水 埋 深 平均 模拟 误差 (MAE ) 
为 0.575 m, 阿 拉 干 断面 上 的 地 下 水 埋 深 平均 模拟 误 
差 (MAE ) 为 0.595 m, 依 干 不 及 麻 断 面 上 的 地 下 水 埋 
深 平 均 模 拟 误 差 (MAE ) 为 0.484 m。 虽 然 这 些 误差 
看 起 来 稍 高 ,但 对 于 河道 日 流量 仅 通 过 年 值 插值 得 
到 的 模型 模拟 来 说 ,模拟 结果 还 是 合理 的 。 另 外 3 个 
断面 上 模拟 与 观测 的 相关 系数 (CC) 也 是 较 高 的 , 它 
们 分 别 是 英 苏 0.783, 阿 拉 干 0.796, 依 干 不 及 麻 
0.945。 总 之 ,通过 与 观测 值 进行 比较 ,进一步 证 实 了 
模型 模拟 的 3 个 断面 的 地 下 水 埋 深 结果 是 合理 的 。 

然后 我 们 分 析 每 个 断面 上 的 平均 地 下 水 埋 深 
的 变化 ,试图 解释 生态 输 水 量 与 断面 地 下 水 埋 深 的 
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输 水 条 件 下 河流 剖 面 地 下 水 埋 深 估算 


变化 模拟 结果 见 图 4。 结 果 表明 ,在 生态 输 水 的 20 a 
里 ,3 个 断面 上 的 地 下 水 位 埋 深 总 体 上 是 逐年 降低 
的 ,对 应 于 地 下 水 水 位 是 逐年 上 升 的 ,这 充分 展示 
了 生态 输 水 对 塔里木 河 下 游 两 岸 的 地 下 水 补给 及 
生态 恢复 起 到 了 积极 的 作用 。 

具体 来 说 ,经 过 20 a 的 生态 输 水 ,塔里木 河 下 游 
的 上 段 的 英 苏 断面 ,中 段 的 阿拉 干 和 下 段 的 依 干 不 
及 麻 的 断面 地 下 水 埋 深 平均 埋 深 从 输 水 前 的 8 mZ 
右上 升 到 4 m 左 右 , 可 见 生 态 输 水 补给 塔里木 河 下 
游 河 岸 地 下 水 的 效果 明显 。 相 比 于 英 苏 断面 上 地 
下 水 位 的 曲折 上 升 ,阿拉 干 和 依 干 不 及 麻 的 地 下 水 
位 是 持续 上 升 的 ,这 主要 与 不 同 断面 上 地 下 水 水 平 
传导 率 的 差异 有 关 。 英 苏 断 面 的 地 下 水 水 平 传导 
率 最 大 ,所 以 地 下 水 埋 深 随 河 道 流量 的 变化 (由 生 
态 输 水 引起 ) 最 明显 。 阿 拉 干 和 依 干 不 及 麻 的 地 下 
水 水 平 传导 率 相 比 英 苏 较 小 ,所 以 在 短 的 输 水 间隔 
期 内 地 下 水 埋 深 受 河道 水 位 的 影响 较 小 ,也 就 是 
说 ,在 输 水 停滞 的 时 间 段 内 ,河道 虽 无 地 表 径 流 , 但 
有 较 浅 的 地 下 水 仍 在 补给 河 两 岸 的 地 下 水 ,最 终 使 
得 地 下 水 位 呈现 为 持续 上 升 的 趋势 。 

比较 3 个 断面 上 地 下 水 埋 深 的 时 间 变 化 可 以 发 
现 , 上 段 英 苏 断 面 的 地 下 水 埋 深 最 浅 , 其 次 是 下 段 
的 依 干 不 及 麻 和 中 段 的 阿拉 干 。 虽 然 中 段 阿 拉 干 
河道 分 配 的 输 水 流量 高 于 下 段 依 干 不 及 麻 ( 约 1.74 


变化 关系 ,3 个 断面 上 1000 m 内 的 平均 地 下 水 埋 深 
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图 3 3 个 断面 上 模拟 与 观测 的 年 均 地 下 水 埋 深 比较 


Fig. 3 Comparisons between simulated and observed mean 


annual groundwater table depths at the three sections 


倍 ) ,但 是 依 干 不 及 麻 的 地 下 水 埋 深 却 低 于 阿拉 十 ， 
该 结果 已 在 刘斌 等 "和 白 元 等 ”的 部 分 年 份 观测 
数据 研究 中 得 到 了 证 明 , 从 侧面 证 明了 拟 二 维 地 下 
水 模型 的 模拟 结果 是 合理 的 。 归 其 原因 与 3 个 断面 
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图 4 3 个 断面 上 模拟 的 年 均 地 下 水 埋 深 变化 
Fig. 4 Variations of simulated mean annual groundwater 


table depths at the three sections 
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上 的 地 下 水 水 平 传导 率 不 同 造成 的 :下 段 依 干 不 及 
麻 的 地 下 水 水 平 传导 率 是 中 段 阿 拉 干 的 2.16 售 。 
高 的 地 下 水 水 平 传导 率 加 快 了 河水 向 河 两 岸 的 地 
下 水 渗透 速度 ,使 得 地 下 水 上 升 的 快 ,最 终 出 现 了 
虽然 中 段 阿 拉 干 断面 上 的 河水 流量 高 于 下 上 段 依 干 
不 及 麻 , 但 下 段 依 干 不 及 麻 的 地 下 水 抬升 高 度 却 高 
于 中 有 段 阿拉 干 的 情况 。 

另外 ,不 同年 份 的 塔里木 河 下 游 生态 输 水 总 量 
有 差异 , 它 是 如 何 影 响 断面 上 地 下 水 埋 深 的 变化 
的 ? 图 5 展示 了 塔里木 河 下 游 生态 输 水 总 量 及 其 3 
个 断面 上 地 下 水 埋 深 的 年 际 变 化 情况 。 由 图 5 可 
知 ,生态 输 水 工程 基本 每 年 都 会 开展 ,并 量 输 水 量 
总 体 呈 增加 的 趋势 ,特别 是 自 2009 年 以 来 年 输 水 量 
增加 明显 。 相 应 地 ,3 个 断面 上 的 地 下 水 埋 深 逐年 
抬升 ,而 且 在 2010 年 以 后 抬升 趋势 明显 。 同 样 地 ， 
地 下 水 开始 下 降 或 放 缓 的 年 份 出 现在 输 水 量 较 少 
的 2007 一 2009 年 和 2014 年 。 可 见 输 水 量 与 地 下 水 
水 位 表现 为 正 相 关 关 系 ,说 明生 态 输 水 对 抬升 下 游 
河岸 的 地 下 水 埋 深 起 到 了 积极 的 作用 。 也 注意 到 ， 
河流 量 减 少 也 会 带 来 河岸 地 下 水 埋 深 的 增 大 ,这 说 
明 塔 里 木 河 下 游 河 岸 要 想 获 得 持续 的 生态 效益 ,也 
tee dee EET EE EY AE ASHIK o 

由 图 5 可 以 看 到 ,地 下 水 埋 深 抬升 (下 降 ) 最 明 


显 的 年 份 ,不 是 对 应 于 输 水 量 最 大 (最 小 ) 的 年 份 ， 
而 是 在 其 之 后 的 1~2 a 内 发 生 。 比 如 英 苏 地 下 水 埋 
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图 5 塔里木 河 下 游 年 生态 输 水 量 及 其 3 个 断面 
上 年 均 地 下 水 埋 深 的 变化 


Fig. 5 Variations of accumulated annual ecological water 


conveyance amounts in the lower reaches of Tarim River and 


mean groundwater table depths at the three sections 


深 局 部 最 大 年 (2005 年 和 2015 年 ) 都 是 在 上 一 年 输 
水 量 出 现 最 低 值 之 后 发 生 的 , 它 在 2018 年 的 最 低地 
下 水 埋 深 是 在 上 一 年 出 现 最 大 输 水 量 之 后 发 生 的 ; 
阿拉 干 和 依 干 不 及 床上 的 地 下 水 埋 深 抬 升 段 
(2004—2006 年 .2012 一 2014 年 .2018 一 2020 年 ) 对 
应 的 年 输 水 量 均 是 下 降 的 ,但 他 们 有 个 共同 的 特点 
是 都 在 前 一 年 出 现 了 最 大 输 水 量 。 由 此 可 见 , 输 水 
量 与 地 下 水 的 年 变化 具有 滞后 性 的 特点 ,这 在 湾 疆 
辉 等 的 和 郭 红 雨 等 加 的 站 点 监测 研究 中 也 被 证 实 ， 
说 明了 拟 二 维 地 下 水 模型 的 地 下 水 变化 模拟 结 
是 合理 的 ,而 且 河 岸 地 下 水 对 生态 输 水 的 响应 具有 
滞后 性 的 特征 。 

除了 比较 不 同 断 面 上 平均 地 下 水 埋 深 的 变化 
外 ,也 比较 了 同一 断面 上 不 同 监测 井 内 地 下 水 埋 深 
变化 。 图 6a~c 分 别 展示 了 英 苏 断面 上 4 个 监测 并 
(C3、C4、C5、C6) 阿拉 干 3 个 监测 井 (G3 .G4、.G5) 和 
依 干 不 及 麻 2 个 监测 井 (H2、H3) 的 年 均 地 下 水 埋 深 
变化 。 结 果 发 现 离 河岸 近 的 监测 井 ( 如 C3、G3 和 
H2) 的 地 下 水 埋 深 经 常 低 于 离 河 岸 远 的 监测 井 ( 如 
C6.G5 和 H3) ,并 且 离 河 近 的 监测 井 地 下 水 埋 深 失 
升幅 度 要 高 于 离 河 远 的 监测 井 。 在 每 个 断面 上 , 监 
测 井 间 的 水 位 差 在 输 水 的 末期 (如 2015 一 2020 年 ) 
相差 不 大 ,这 充分 说 明 多 年 河道 生态 输 水 对 河岸 地 
下 水 埋 深 的 抬升 起 到 了 积极 作用 。 这 些 结果 说 明 
拟 二 维 地 下 水 模型 模拟 同一 个 断面 上 不 同 监 测 井 
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图 6 断面 监测 井 的 年 均 地 下 水 埋 深 变化 
Fig. 6 Variations of mean annual groundwater table 


depths at different monitoring wells of each section 
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图 7 3 个 断面 上 模拟 的 地 下 2 m 深 平均 
土壤 湿度 的 年 际 变化 


Fig. 7 Annual variations of simulated mean 2-meter-depth 


soil moistures at the three sections 


的 地 下 水 埋 深 变化 是 合理 的 。 
2.3 断面 土壤 湿度 的 长 期 模拟 

塔里木 河 下 游 生态 输 水 目的 不 仅 是 抬升 地 下 
水 埋 深 ,更 重要 的 是 提高 表层 土壤 中 非 饱 和 土壤 水 
含量 (也 称 土壤 湿度 ) ,满足 河岸 植被 的 生长 需求 和 
改 恨 土壤 盐 渍 化 的 条 件 。 为 此 , 非 愧 和 土壤 水 含量 
在 输 水 条 件 下 的 变化 也 是 本 研究 关注 的 一 个 重 
点 。 由 于 3 个 断面 的 地 表 处 于 不 同 的 海拔 高 度 ,为 
了 统一 起 见 ,我 们 选取 地 下 埋 深 2m 内 的 平均 土壤 
湿度 模拟 结果 进行 分 析 比 较 , 最 终 20 a 的 土壤 湿度 
模拟 结果 如 图 7 所 示 。3 个 断面 上 的 土壤 湿度 在 20 a 
的 输 水 过 程 中 是 逐年 增加 的 ,其 值 分 别 由 最 初 的 
0.20 上 升 到 0.30 以 上 ,其 中 英 苏 和 依 干 不 及 麻 的 上 
升值 都 在 0.35 左右 。3 个 断面 上 ,上 上段 英 苏 的 20 a 
平均 土壤 湿度 最 高 (均值 为 0.294) ,其 次 是 下 段 依 干 
不 及 麻 (0.29) 和 中 段 阿拉 干 (0.26)。 另 外 ,土壤 湿 
度 变 化 起 伏 最 大 的 英 苏 断面 ,其 它 两 段 面 上 升 比较 
平稳 。 这 些 变化 特点 都 与 各 自 的 地 下 水 埋 深 变化 
趋势 比较 相似 ,可 以 算出 20a 内 地 下 水 埋 深 变化 与 
土壤 湿度 变化 的 相关 系数 分 别 是 英 苏 0.892 ,阿拉 干 
0.867 和 依 干 不 及 麻 0.957。 由 于 垂直 的 土壤 饱和 传 
导 率 值 在 3 个 断面 上 是 相等 的 ,所 以 地 下 水 与 土壤 
水 间 呈 高 相关 关系 ,可 以 认为 土壤 湿度 的 增加 主要 
是 由 地 下 水 埋 深 抬升 所 引起 的 。 

土壤 湿度 随 生 态 输 水 量 的 变化 结果 见 图 8。 正 
如 土壤 湿度 的 变化 和 地 下 水 变化 具有 和 较 高 相关 性 
一 样 , 土 壤 湿 度 与 生态 输 水 量 的 关系 也 满足 2 个 特 
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图 8 塔里木 河 下 游 生态 输 水 量 与 其 3 个 断面 上 
2 m 深 平均 土壤 湿度 的 年 变化 


Fig. 8 Annual variations of accumulated ecological water 


conveyance amounts in the lower reaches of Tarim River and 


mean 2-meter-depth soil moistures at the three sections 


点 : 正 相 关 性 和 滞后 性 。 土 壤 湿 度 随 着 输 水 时 间 和 
输 水 量 的 增加 而 增加 ,特别 是 后 10 a, 随 着 输 水 量 的 
大 幅 增加 ,土壤 湿度 相 比 前 10a 有 明显 的 增加 。 土 
壤 湿 度 与 生态 输 水 量 的 滞后 性 也 非常 明显 ,比如 3 
个 断面 在 生态 输 水 量 最 高 年 (2003、2011、2017 年 ) 
之 后 的 2~3 a 内 ,土壤 湿度 仍然 表现 为 增加 的 趋势 ; 
在 生态 输 水 最 低 年 (2004、2008、2014 年 ) 对 应 的 土 
襄 湿 度 都 不 是 局 部 最 小 值 ,而 是 在 其 后 1~2 a 内 出 
现 了 局 部 最 小 值 。 


3 结论 


本 文 利用 发 展 的 包含 土壤 水 和 地 下 水 运动 的 
拟 二 维 地 下 水 模型 对 塔里木 河 下 游 输 水 条 件 下 上 、 
中 、 下 段 的 3 个 断面 ( 英 苏 、 阿 拉 干 和 依 干 不 及 麻 ) 地 
下 水 埋 深 和 土壤 湿度 进行 了 20 a 长 期 模拟 。 通 过 率 
定期 和 11 a(2010 一 2020 年 ) 的 地 下 水 埋 深 模 拟 结 
与 观测 数据 比较 ,发 现 两 者 是 较 一 致 的 ,表明 模型 
的 模拟 结果 是 合理 的 ,从 而 证 明了 该 拟 二 维 地 下 水 
模型 在 塔里木 河 下 游 断面 地 下 水 模拟 上 的 适用 
性 。 然 后 ,根据 合理 模拟 的 3 个 断面 20a 地 下 水 埋 
深 的 结果 ,横向 分 析 每 个 断面 上 在 生态 输 水 条 件 下 
的 20a 地 下 水 埋 深 变化 ,纵向 分 析 3 个 断面 上 地 下 
水 埋 深 变 化 及 其 与 生态 输 水 量 的 关系 。 另 外 ,对 断 
面 土壤 湿度 的 相应 变化 也 进行 了 分 析 。 主 要 结论 
如 下 : 

(1) 为 了 提高 模型 在 各 个 断面 上 对 地 下 水 动态 
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模拟 的 适用 性 ,首先 对 拟 二 维 地 下 水 模型 在 断面 上 
的 模拟 进行 参数 率 定 。3 个 断面 上 的 模型 率 定 参数 
均 为 地 下 水 水 平 传导 率 。 英 苏 断 面 模拟 的 率 定期 
为 第 2 次 生态 输 水 ,阿拉 干 和 依 干 不 及 麻 断 面 模拟 
的 率 定期 为 第 4 次 生态 输 水 。 这 样 的 设置 主要 是 根 
据 断 面 地 下 水 观测 数据 的 密集 程度 而 确定 。 最 终 
英 苏 、 阿 拉 干 和 依 干 不 及 麻 断 面 上 模型 率 定 的 地 下 
水 最 优 模拟 的 平均 绝对 值 误差 (MAE) 分 别 为 
0.194 m .0.214 m、0.189 m 并 且 每 个 断面 上 监测 井 地 
下 水 位 离 河 岸 越 远 而 越 低 ,说 明了 拟 二 维 地 下 水 模 
型 对 塔里木 河 下 游 输 水 条 件 下 河岸 断面 地 下 水 位 
模拟 是 合理 的 ,可 以 在 下 游 河 岸 进行 长 期 模拟 。 
(2) 对 最 优 的 拟 二 维 地 下 水 模型 ,分 别 在 3 个 
斯 面 上 进行 20 a 输 水 条 件 下 的 地 下 水 埋 深 的 模拟 ， 
用 2010 一 2020 年 的 地 下 水 埋 深 数据 验证 模拟 结 
结果 表明 每 个 断面 的 地 下 水 埋 深 年 均 模 拟 误 差 
(MAE) 2) 460.5 m 左右 ,这 在 缺乏 准确 的 日 河水 流 
量 数据 情况 下 ,模拟 结果 被 认为 是 合理 的 ,从 而 证 
明了 模型 对 塔里木 河 下游 断 面 地 下 水 长 期 的 模拟 
是 合理 的 。 然 后 结合 20 a 的 模拟 结果 进行 分 析 , 表 
明 随 着 输 水 年 份 的 增加 ,3 个 断面 上 的 地 下 水 埋 深 
均 有 明显 的 抬升 ,地 下 水 位 埋 深 从 输 水 前 的 8 m 左 
右 抬升 到 输 水 后 的 4m 左 右 。 特 别 是 2009 年 之 后 ， 
随 着 年 输 水 量 的 增加 ,断面 上 的 地 下 水 埋 深 抬升 明 
显 ,这 些 将 会 对 于 河岸 的 生态 植被 恢复 起 到 了 积极 
的 作用 。 下 段 依 干 不 及 麻 的 地 下 水 传导 率 参 数值 
是 中 段 阿 拉 干 的 2.16 倍 ,这 使 得 虽然 中 段 河道 流量 
高 于 下 有 段 河道 流量 ,但 最 终结 果 是 下 有 段 依 干 不 及 麻 
的 地 下 水 水 位 高 于 中 有 段 阿拉 干 ,其 中 地 下 水 水 平 传 
导 率 参数 起 着 重要 的 作用 。 通 过 比较 塔里木 河 下 
游 输 水 量 与 3 个 断面 上 地 下 水 埋 深 的 变化 关系 ,发 
现 同一 断面 上 地 下 水 埋 深 对 河水 流量 有 一 定 的 沾 
后 正 相 关 作 用 。 另 外 ,也 分 析 了 同一 断面 不 同 监测 
井 的 地 下 水 埋 深 模拟 变化 ,发现 离 河 近 的 监测 井 地 
下 水 埋 深 抬升 幅度 高 于 离 河 远 的 监测 井 。 这 些 结 
果 从 纵向 (不 同 断 面 ) 和 横向 (同一 断面 ) 上 ,证实 了 
模型 对 输 水 条 件 下 河岸 地 下 水 长 期 模拟 的 合理 性 。 
同时 , 输 水 量 与 地 下 水 水 位 的 正 相关 与 沾 后 性 ,说 
明了 进行 生态 输 水 对 抬升 下 游 河 岸 的 地 下 水 埋 深 
起 到 了 积极 的 作用 ,并 且 塔 里 木 河 下 游 河岸 要 想 获 
得 持续 的 生态 效益 ,需要 提供 间 欣 性 的 生态 输 水 。 


(3) 在 缺乏 观测 的 情况 下 ,也 对 模拟 断面 2 m 
深 的 土壤 湿度 进行 了 分 析 。 结 果 发 现 , 随 着 生态 输 
水 量 的 增加 ,土壤 湿度 也 有 明显 的 增加 ,其 值 从 输 
水 前 的 0.20 增 加 到 20 a 输 水 后 的 0.30 以 上 ,其 中 英 
苏 和 依 干 不 及 麻 的 上 升值 都 在 0.35 左右 ,这 些 对 于 
河岸 的 生态 植被 恢复 起 到 了 积极 的 作用 。 土 壤 水 
与 地 下 水 位 的 变化 趋势 有 较 好 的 正 相关 性 ,两 者 相 
关系 数 在 0.85 以 上 。 由 于 模型 在 垂直 土壤 水 传导 
率 参 数 上 取 值 一 样 , 因 此 可 以 推断 土壤 湿度 的 变化 
主要 是 由 地 下 水 埋 深 抬 升 引 起 。 通 过 分 析 输 水 量 
与 土壤 湿度 的 年 变化 ,发 现 土壤 湿度 与 输 水 量 间 也 
有 滞后 性 特征 。 

本 研究 在 用 模型 解决 输 水 条 件 下 地 下 水 动态 
变化 实际 问题 上 ,用 参数 率 定 法 消除 了 模型 在 描述 
实际 问题 的 误差 ,但 这 毕竟 是 有 限 的 。 如 果 能 把 更 
多 的 物理 过 程 如 植被 光合 作用 土壤 的 冻 融 和 不 同 
植被 动态 变化 等 加 入 到 模型 描述 中 ,将 会 得 到 好 的 
模型 模拟 结果 。 另 外 还 需要 引入 多 种 数据 源 来 对 
模型 做 进一步 的 发 展 与 验证 ,比如 收集 实测 的 断面 
河水 流量 或 水 位 数据 .地表 降 水 .实际 蒸 散 发 量 和 
更 多 的 监测 井 地 下 水 埋 深 等 数据 ,对 当地 的 土壤 质 
地 参数 做 多 样本 点 的 实测 ,收集 更 准确 的 土壤 湿度 
数据 等 ,只 有 这 样 才能 获得 更 加 准确 的 地 下 水 和 土 
壤 水 模拟 ,进而 建立 准确 的 河水 流量 与 地 下 水 位 动 
态 关 系 式 。 最 终 根据 河岸 植被 生态 所 需 水 量 ,设计 
不 同情 景 的 输 水 策略 ,比如 选择 不 同 的 输 水 期 和 输 
水 量 等 ,能 以 最 小 的 输 水 量 最 大 程度 地 满足 塔里木 
河 下 游 的 生态 环境 安全 ,这 对 于 干旱 区 生态 水 资源 
管理 具有 重要 的 意义 和 参考 价值 。 
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Estimation of riparian groundwater table depth in the lower reaches of Tarim 


River under long-term water conveyance 
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Abstract: Accurate estimation of the groundwater table dynamic variation in a riparian zone under ecological wa- 
ter conveyance conditions is of great importance for the sustainable development of water resources management 
in arid areas. Based on the data of river discharges from ecological water conveyance and the groundwater tables 
of riparian monitoring wells from 2001 to 2020, this study simulated the dynamic variation of groundwater table 
at three sections of Yingsu, Alagan, and Yiganbujima in the lower reaches of the Tarim River, Xinjiang, China us- 
ing a developed quasi- two-dimensional groundwater model, including the groundwater and soil water motion 
equations. In the calibration period, the observed groundwater table data was used to adjust the critical parameter 
(1.e., horizontal hydraulic conductivity) of the quasi-two-dimensional model, and the average simulation error of 
the monitoring wells at each of the three sections was approximately 0.2 m. In the validation period, the ground- 
water table data from 2011 to 2020 at the monitoring wells of the three sections were compared with the simulat- 
ed results, and the average error at each of the three sections was approximately 0.5 m. Therefore, the simulated 
results of the groundwater table were acceptable under the lack of accurate river level data and the model was 
suitable for simulating the variation of the riparian groundwater table under ecological water conveyance in the 
lower reaches of the Tarim River. Thus, based on the simulation results from 2001 to 2020, the groundwater table 
and soil moisture variations were analyzed to demonstrate the response of riparian groundwater and soil moisture 
to the ecological water conveyance project. The results show that the groundwater table and soil moisture at the 
three sections evidently increased in 20 years of water conveyance periods. Specifically, the groundwater table in- 
creased from approximately 8 m before ecological water conveyance to nearly 4 m, and the corresponding soil 
moisture rose from 0.20 to 0.35 above. The increasing trend is more significant with the substantial increase of 
the annual water conveyance amount since 2009. Moreover, the variation of groundwater table presented a certain 
lag following the change in the water conveyance amount, and a positive correlation between soil moisture and 
groundwater table makes the soil moisture lags behind the water conveyance amount. Compared with the river 
discharge, the groundwater hydraulic conductivity plays a more important role in the variation of groundwater ta- 
ble at cross-section; thus, increasing the upswing amplitude of the groundwater table at Yiganbujima sections of 
the lower segment compared with that of the Alagan section of the middle segment in lower reaches of the Tarim 
River. Additionally, the inter- annual variations of groundwater and water conveyance amount also indicate that 
the intermittent ecological water conveyance is essential to achieve a sustainable ecological benefit in the lower 
reaches of the Tarim River. 

Key words: quasi-two-dimensional groundwater model; ecological water conveyance; long-term water table 


simulation; the lower reaches of Tarim River 


